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Abstract. We study the influence of hadron in- 
teraction features on the high-energy atmospheric 
neutrino spectrum. The ID calculation is performed 
with use of the known high-energy hadronic models, 
SIBYLL 2.1, QGSJET-II, Kimel and Mokhov, for 
the parameterizations of primary cosmic ray spectra 
issued from the measurement data. The results are 
compared with the Frejus data and AMANDA-II 
measurements as well as with other calculations. 
Sizable difference of the neutrino fluxes (up to the 
factor of 1.8 at 1 TeV) that are obtained with 
the SIBYLL and QGSJET-II appears to be rather 
unexpected keeping in mind the hadron and muon 
flux calculations in the same energy region. 
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I. Introduction 

Atmospheric neutrinos (AN) appear in decays of 
mesons (charged pions, kaons etc.) produced through 
collisions of high-energy cosmic rays with air nuclei. 
The AN flux in the wide energy range remains the 
issue of the great interest since the low energy AN 
flux is a research matter in the neutrino oscillations 
studies, and the high energy atmospheric neutrino flux 
is now appearing as the unavoidable background for 
astrophysical neutrino experiments 0, @, 0, 0, 0, 
0. To present day a lot of calculations are made of 
atmospheric neutrino fluxes, among which 0, 0, 0, 

una, ed, ma, in, m, ma (see ai so ma, ed, ma, 

[19], [20 1 for a review of ID and 3D calculations of the 
AN flux in the wide energy range). 

In this work we present results of new one- 
dimensional calculation of the atmospheric muon neu- 
trino flux in the range 10-10 7 GeV made with use of the 
hadronic models QGSJET-II 03 ED, SIBYLL 2.1 ||22) 
as well as Kimel and Mokhov (KM) [23 ] that were tested 
also in recent atmospheric muon flux calculations 11241 . 
We make an attempt to learn how strongly the diversity 
of hadronic interaction models influences on the high- 
energy spectrum of atmospheric neutrinos. 

II. The method and input data 

The calculation is performed on the basis of the 
method 11251 of solution of the hadronic cascade equa- 
tions in the atmosphere, which takes into account non- 
scaling behavior of inclusive particle production cross- 
sections, the rise of total inelastic hadron-nuclei cross- 



sections, and the non-power law primary spectrum (see 
also Il24l ). As the primary cosmic ray spectra and 
composition in wide energy range used is the model 
recently proposed by Zatsepin and Sokolskaya (ZS) ll26l . 
which fits well the ATIC-2 experiment data 11271 and sup- 
posedly to be valid up to 100 PeV. The ZS proton spec- 
trum at E > 10 6 GeV is compatible with KASCADE 
data 11281 as well the helium one is within the range of 
the KASCADE spectrum recontructed with the help of 
QGSJET 01 and SIBYLL models. Alternatively in the 
energy range 1 — 10 6 GeV we use the parameterization 
by Gaisser, Honda, Lipari and Stanev (GH) Ifl8l, 11291, 
the version with the high fit to the helium data. Note 
this version is consistent with the data of the KASCADE 
experiment at Eo > 10 6 GeV that was obtained (through 
the EAS simulations) with the SIBYLL 2.1. 

To illustrate the distinction of the hadron models em- 
ployed in the computations, it is appropriate to compare 
the spectrum- weighted moments (Table computed for 
proton-air interactions (for 7 = 1.7): 



(E ) = 



x 1 da v 



'pA 



dx 



dx, 



(1) 



where x = E c /Eo, c — p, n 1 7r ± , AT ± . The values in 
Table [Qdisplay approximate scaling law both in SIBYLL 
2.1 and KM and little violation of the scaling in the 
QGSJET-II for p and tt ± . 

III. Atmospheric muon neutrino fluxes 

Along with major sources of the muon neutrinos, 71^2 
and K^2 decays, we consider three-particle semileptonic 
, the contribution originated from de- 



decays, K^ 3 , 
cay chains K 



— > 7T 



(K° s 



K 



T^TT ) 



as well as small fraction from the muon decays. 

One can neglect the 3D effects in calculations of the 
atmospheric muon neutrino flux near vertical at energies 
E > 1 GeV and at E > 5 GeV in case of directions 
close to horizontal (see [19], l20l ). Fractions of the 
neutrino flux near vertical from pion and kaon decays 
are shown in Fig. Q] These calculations are made for the 
model primary spectrum by GH [18] (Fig. [Taji as well for 
the model by ZS [26] that comprises the results of ATIC- 
2 experiment 11271 (Fig. [Tbl . Note the similar ratio for 
muon fluxes differs from that of neutrino fluxes: at 10 3 
GeV the [IkI ^ ratio is about 0.25, while the vk /v^ is 
about 4 (see also Fig. 4 in Ref. Ifl8l ). 
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TABLE I: Spectrum weighted moments z pc (Eo), calculated for 7 = 1.7 



Model 


Eo, GeV 


Zpp 


2pn 


Z p-K + 


z pn~ 


z pK+ 


z pK~ 




10 a 


0.174 


0.088 


0.043 


0.035 


0.0036 


0.0030 


QGSJET-II 


10 3 


0.198 


0.094 


0.036 


0.029 


0.0036 


0.0028 




10 4 


0.205 


0.090 


0.033 


0.028 


0.0034 


0.0027 




10 2 


0.211 


0.059 


0.036 


0.026 


0.0134 


0.0014 


SIBYLL 2.1 


10 3 


0.209 


0.045 


0.038 


0.029 


0.0120 


0.0022 




10 4 


0.203 


0.043 


0.037 


0.029 


0.0097 


0.0026 




10 M 


0.178 


0.060 


0.044 


0.027 


0.0051 


0.0015 


KM 


10 3 


0.190 


0.060 


0.046 


0.028 


0.0052 


0.0015 




10 4 


0.182 


0.052 


0.046 


0.029 


0.0052 


0.0015 
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(a) Calculation for the GH primary spectrum 1181 
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(b) Calculation for Zatsepin and Sokolskya model |26| 



Fig. 1: Fraction of {y^ 



flux from kaons (solid lines) and pions (dashed) calculated for 6 = 0° 
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Fig. 2: Ratio of the and fluxes calculated with 
KM model for GH and ZS primary spectra. 

The ratio v^jv^ calculated with KM model for the 
two primary spectra, GH and ZS, is plotted in Fig. [2] 
The wavy shape of the ratios apparently visible in 
Figs. [TBI andl2"lreflects the peculiarities of the ZS spectra. 

A comparison of (y^ + v^) flux calculations for 
the three hadronic models under study is made 
in Table [TTJ column 1, 2 and 3 presents the 
flux ratio, 0™^^, <^ GSJET - n) Mr } ^d 
</>^ IBYLL ' > /</>i^ GSJET n ^ correspondingly, calculated at 



TABLE II: Ratio of the fluxes at 8 = 0° (90°) 
calculated with the SIBYLL 2.1, QGSJET-II, and KM 



E v , GeV 


1 


2 


3 






GH 




10 2 


1.65 (1.22) 


0.97 (0.85) 


1.70 (1.44) 


10 3 


1.71 (1.46) 


0.96 (0.92) 


1.78 (1.59) 


10 4 


1.60 (1.57) 


0.96 (0.96) 


1.67 (1.64) 


10 5 


1.54 (1.49) 


0.99 (0.96) 


1.56 (1.55) 






ZS 




10 2 


1.58 (1.26) 


1.00 (0.91) 


1.58 (1.38) 


10 3 


1.64 (1.39) 


0.95 (0.92) 


1.73 (1.51) 


10 4 


1.55 (1.46) 


0.96 (0.95) 


1.61 (1.54) 


10 5 


1.37 (1.23) 


0.91 (0.83) 


1.51 (1.48) 


10 6 


1.10 (0.95) 


0.61 (0.55) 


1.80 (1.73) 


10 7 


0.89 (0.75) 


0.48 (0.43) 


1.85 (1.74) 



6 = 0° and 90° (in brackets) with usage of the GH 
and ZS primary spectrum. One can see that usage 
of QGSJET-II and SIBYLL models leads to apparent 
difference of the muon neutrino flux, as well as in the 
case of SIBYLL as compared to KM (unlike the muon 
flux, where SIBYLL and KM lead to very similar results 
ED). On the contrary, the QGSJET-II neutrino flux is 
very close to the KM one: up to 100 TeV the difference 
does not exceed 5% for the GH spectrum and 10% for 
the ZS one at 6 = 0°. While the muon flux discrepancy 
in the QGSJET-II and KM predictions is about 30% at 
vertical ll24ll . The origin of differences is evident: the 
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Fig. 3: Comparison of the two independent calculations 
for the GH spectrum. 

kaon production ambiguity. 

Figure [3] shows this work calculations of the neu- 
trino flux (lines) in comparison with the result of Barr, 
Gaisser, Lipari, Robbins and Stanev (BGLRS) fl]0) ob- 
tained with use of the TARGET 2.1 (symbols). All these 
computations are performed for the GH primary spectra. 
One can see the calculations for KM and TARGET 2.1 
are in close agreement in the range 10 — 10 4 GeV (near 
horizontal) as well as at E v < 200 GeV near vertical. 

In Fig. [4] presented is the comparison of different 
calculations of the AN flux along with the data of the 
AMANDA-II EL El and Frejus (30) experiments. More 
comparisons of the low and high energy flux calculations 
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Fig. 4: Comparison of the present calculation as well as 
the previous ones (by Volkova [7 1, Butkevich et al. IflOl . 
Lipari HD, Naumov et al. ||T3|) with the data of the 
AMANDA-II [4| and Frejus ESI experiments. This work 
calculation codes are in the right top corner. 



may be found in Refs. O, ED, ED, (T7J, Ifl9l . 

Figure [5] presents the comparison of this work calcu- 
lation of the conventional (from fi, n, i^-decays) and 
prompt muon neutrino flux with some of previous 
ones Q3), Q7), ED, (32), (33]. The conventional flux 
here was calculated with use of QGSJET-II model 
combined with the Zatsepin and Sokolskaya primary 
spectrum (thin lines). Dashed lines mark the calculation 
by Naumov, Sinegovskaya and Sinegovsky ITPJI . IfTTll of 
the conventional muon neutrino fluxes for 9 = 0° and 
90°. Bold dotted line (curve 1) shows the sum of the 
prompt neutrino flux by Volkova and Zatsepin (VZ) l32l 
and the conventional one due-to the QGSJET-II + ZS 
model at 9 = 90°. Dash-dotted line (2) marks the sum 
of the QGSJET-II conventional flux (9 = 90°) and the 
prompt neutrino contribution due to the recombination 
quark-parton model (RQPM) ll3D . Solid line 4 shows the 
same for the prompt neutrino flux due to the quark-gluon 
string model (QGSM) ED (see also H3, E), ll34l ). 
Also shown are the two of the prompt neutrino predic- 
tions by Gelmini, Gondolo and Varieschi (GGV) 1331 
(curves 3 and 5). 

Notice that recent evaluation of the prompt neutrino 
flux obtained with the dipole model (DM) E3, is rather 
close to the QGSM prediction at E v > 10 6 GeV, 
keeping in mind that the theoretical uncertainty absorbs 
a difference of the DM and QGSM fluxes. 

The prompt neutrino fluxes at E v = 100 TeV are 
presented in Table [III] along with the upper limit on 
the astrophysical muon neutrino diffuse flux obtained in 
AMANDA-II experiment [4|. Note that the QGSJET- 
II+GH flux appears to be the most low flux of the 
conventional atmosperic neutrinos at high energies. 

TABLE III: Atmospheric neutrino flux at E v = 100 TeV 
vs. the AMANDA-II restriction for the + flux 



Model 


E'£(t> u , (cm 2 ssrr 1 GeV 


conventional + P M : 


0° 90° 


QGSJET-II + ZS 


1.20 x 10~ 8 10.5 x IO" 8 


QGSJET-II + GH 


1.11 x IO" 8 9.89 x 10~ 8 


prompt Ufj, + P M : 


90° 


vz 03 


8.12 x IO" 8 


RQPM |31 1 


4.61 x IO" 8 


QGSM 1 3 1 1 


1.22 x IO" 8 


AMANDA-II upper limit |4| 


7.4 x 10" 8 



IV. Summary 

The calculations of the high-energy atmospheric muon 
neutrino flux demonstrate rather weak dependence on 
the primary specrtum models in the energy range 10 — 
10 5 GeV. However the picture seems to be less steady 
because of sizable flux differences originated from the 
models of high-energy hadronic interactions. As it can 
be seen by the example QGSJET-II and SIBYLL 2.1, the 
major factor of the discrepancy is the kaon production 
in nucleon-nucleus collisions. 

A common hope that atmospheric muon fluxes might 
be reliable tool to promote the discrimination between 
the hadronic interaction models seems to be rather 
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Fig. 5: Fluxes of the conventional and prompt muon neutrinos along with data points from the AMANDA-II (4) 
and Frejus [ 30 1 experiments. Codes of curves marking the prompt neutrino flux at 8 — 90° are folows: 1 - VZ [32|; 
2 - RQPM OP ; 3 - GGV 1331 (the case of A = 0.5, where A is exponent of the gluon distribution at low Bjorken 
x) ; 4 - QGSM [31 1; 5 - GGV (A = 0.1). Curves just below the 3, 4 and 5 display the coresponding flux at = 0°. 



illusive as the key differences in the ir/K production 
impact variously on the neutrino flux and muon one. For 
the high-energy neutrino production at the atmosphere 
the kaon yield in nucleon-nucleus interactions is more 
strong factor in comparison with that for production 
of the atmospheric muons, despite on their common to 
neutrinos origin. 

Inasmuch as the atmospheric prompt neutrino flux 
weakly depends on the zenith angle (near 100 TeV), one 
may refer the AMANDA-II restriction just to the prompt 
neutrino flux model. Thus one may consider both RQPM 
and QGSM to be consistent with the AMANDA-II upper 
limit for diffuse neutrino flux. 
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